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ABSTRACT 


A nuclear electric propulsion concept using a thermionic reactor 
inductively coupled to a magnetoplasmadynamic accelerator (MPD tire jet) 
is described and the results of preliminary analyses are presented. In 
this system, the MPD thrvister operates intermittently at hi^er voltages 
and power levels than the thermionic generating vinit. A typical thrust 
pulse from the MPD arc jet in this study is characterized by power levels 
of 1 to MWe, a duration of 1 msec, and a duty cycle of '^ 20 %. The 
thermionic generating unit operates continuously but with a lower power 
level ~0.4 MWe. Energy storage between thrust pulses is provided by 
building up a large current in an inductor using the output of the 
thermionic converter array. Periodically, the charging current is 
interrupted and the energy stored in the magnetic field of the inductor 
is utilized for a short duration thrust pulse. The results of the pre- 
liminary analysis show that a coupliog effectiveness of approximately 
85 to 90 % is feasible for a nominal UOO KWe system with an inductive 
unit suitable for a flight vehicle. Optimized values of the total 
specific mass of the system including the thermionic reactor, the in- 
ductor, and the MPD thruster are estimated in the range of 23 to 2U kg/KWe. 
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1.0 NOMENCLATURE 


Roman Letters 

B : Parameter used in calculating the inductance of a long solenoid 

s 

d : Mean diameter of the inductor, cm 

F : Parameter used in calculating the inductance of a long solenoid 

I : Cwrent, amps or kiloamps (kA) 

I : Maximum value of the current in a cycle, amps or kiloan^js (kA) 

I : Minimum value of the current in a cycle, amps or kiloamps (kA) 

I : Short circuit c\u*rent of thermionic converter array, amps 

o 

L : Inductance of the solenoid inductor, ph 

h : Length of the inductor, cm 

Mj^ ; Mass of the inductor, kg 
N : Number of turns in the inductor 

P : Net power supplied to thruster, KWe 

P : Maximum electric power available from the thermionic reactor, KWe 

o 

R : Resistance of the inductor winding, ohms 

R : Total resistance of the current buildup part of the circuit, ohms 

P 

R : Resistance of the MPD thruster part of the circuit, ohms 

s 

s : Thickness of the inductor winding, cm 

T : Time between pulses for current buildup, msec 

t* : Dxiration of the thrust pulse, msec 

V : Voltage of the thermionic converter array, volts 

c 

V : Parameter in the approximate thruster V-I curve, 30 volts 

F 

V : Open circuit voltage of the thermionic converter array, volts 

o 

: Voltage of the MPD thruster, volts 


Greek. Leters: 


: Specific mass of the MPD thruster, kg/KWe 

a : Specific mass of the complete system, kg/KWe 

O 

a_ : Specific mass of the thermionic reactor, kg/KWe 

R 

: Specific mass of the switching transistors, kg/KWe 

^TR ’ Specific mass of the transistor radiator, kg/KWe 

$ : Parameter in the approximate thruster V-I curve, 1.2 

n : Coupling effectiveness, % 

T : Time constant for current buildup, msec 

p : Temperature dependent resistivity, yJ^-cm 

9 : Temperatvire of the inductor, °K 


2.0 INTRODUCTION 


Providing propulsion power for deep space missions is often best accom- 
plished with a nuclear electric system. An attractive choice for the power 
system is a nuclear reactor which uses thermionic conversion for electrical 
generation. The thermionic reactor is a compact, simple system with no moving 
parts and a hi^ heat rejection temperature which permits a smaller radih,"^or to 
be used for dissipating the waste heat. Other advantages for spacecraft power 
include reliability (no single point failures), long life, and a low specific 
mass (Ref. l). 

The characteristics of the thruster tinit must be optimized fpr each mission. 
The thrust variables which must be considered include efficiency of conversion, 
specific impulse, power level, and specific mass. The MPD thruster has many 
advantages for spacecraft electric propulsion. A principal advantage of MPD 
devices over ion rockets is the high thrust density available with the arc jet 
(typically 1 to 10 N/cm ), This high thrust density usually allows a signifi- 
cant savings in the mass of the system. 

However, efficient operation of the MPD thruster is achieved at power levels 
(in the megawatt range) which are hi^er than optimum for many missions. Further- 
more, the output of the thermionic converter array is low voltage DC which re- 
quires voltage transformation to match the reqviirements of the MPD arc jet. The 
differences in power levels between the thruster and the thermionic generating 
vinit combined with the need for power conditioning have led to a system design 
in which power is transferred to the thruster by means of inductive energy storage. 

The thruster is operated in quasi-steady or pulse mode while the thermionic 
generating unit operates continuously. Between pulses, energy is stored in the 
magnetic field of a large inductor and dissipeted in each short duration thrust 
pulse. The operating cycle begins with a relatively long charging period (sev- 
eral msec). During this time the current is built up in the windings of the 
inductor. The charging current is interrupted at the end of the charging period. 
Next follows the short thrust period (~ 1 msec) where the current decays rapidly 
to a lower value. The rapidly changing magnetic flux inihe inductor generates 
a high voltage to power the thruster. The only other required power conditioning 
element is the switching device(s) to interrupt the charging current- 


3.0 SYSTEM DESCRIPTION 

3.1 Vehicle Configuration 

The thermionic reactor design and spacecraft config\iration which 
form a has is of this work Lave heen developed by previous studies at Jet 
Propulsion Laboratory (Refs. 1 & 2). The design chosen for the spacecraft 
is an end thrust type of vehicle which is shown in Fig. 1. The thermionic 
reactor is located at the rear of the spacecraft behind the radiation shield. 
The inductor must be configured to fit inside the 23° shadow cone of the 
radiation shield. 

3.2 Thermionic Reactor 

Design of the thermionic reactor has been carried out primarily by 
personnel at Los Alamos Scientific Laboratory as a subcontract to JPL (Refs. 

2 & 3). The reactor system is an out-of-core design which uses high tenper- 
ature heat pipes to transfer heat from a fast spectr'am compact reactor to an 
array of thermionic converters. The configuration of the reactor system 
(Ref. 2) is shown in Fig. 2. Electrical connections between converters are 
made in a direction transverse to the axis of the heat pipes. The output bus 
bars are arranged in 6 layers which are accessible for connections at the peri- 
phery of the thermionic converter array. 

The emitter temperature of the thermionic converters is l650°K and 

the collector temperature is ~ 900°K. Average power density is approximately 

2 2 
6 watts/cm . The area of each converter is a l62.6 cm . There are 5^0 con- 
verters in the complete power package (Ref. 2). 

Various combinations of series-parallel connections for the therm- 
ionic converter array can produce output voltages from ~ 3 volts to ~ 5^ volts 
with correspondingly different currents. The optimum voltage for the inductive 
coupled MPD thruster depends upon resistive and switching losses as well as the 
variation of thruster performance with current. Determination of the optimum 
connection arrangement is not possible because detailed data on the variation 
of thruster efficiency with operating cxirrent are not available. Compromise 
values of the reactor current and voltage were selected for this study which 
produce currents of 12,000 amps at ~ 39 volts. The total power at the bus 
bars is Kl/e. 
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Fig. 2 JPL/LASL Out-of-Core Thermionic Reactor Design 
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Tie current- .Itage characteristic of a typical thermionic device 
is shown in Fig. 3a* For purposes of the analysis, a linear approximation 
was used for the cvtrrent-voltage characteristic of the entire thermionic array 
as is shown in Fig. 3b. The open circuit /oltage, shown on the figure as V^, 
is 78 volts and the short circuit current, I^, is 2l*,390 an^s. The approximate 
I-V characteristic is written as 

= V^(l - I/I ). (1) 

coo 

The maximuffi power from the reactor, occurs at approximately V^/2 and 1^/2; 

i.e. 


P 

o 


I V 
o o 

k 


(2) 


3 . 3 Magnetoplasmadynamic Thruster 

Research on the quasi-steady MPD thruster is being carried out by 
the School of Engineering and Applied Sciences at Princeton University (Ref. U). 
This Princeton program includes experiments to evalioate the thruster character- 
istics and demonstrate the use of inductive coupling. In this set of experi- 
ments, the power source is simulated by a capacitor bank. 

The MPD thruster consists of a two electrode device and a propellant 
injection system as shown schematically in Fig. h. The cylindrical cathode at 
the center is surrovinded by an annular annode. Propellant is injected throu^ 
the insulating beck plate and accelerated out as a plasma to produce thrust. 
Current traveling axially down the cathode produces an azimuthal magnetic field. 
When the current then flows through the plasma it produces a 7 x B acceleration 
to move the plasma out the back of the thruster. 

High currents are required by the thruster to produce the aximuthal 
self magnetic field. Operation of an MPD device at lower currents is possible 
but usually requires an externally generated magnetic field. The self field 
tj'pe operation is much preferred for better efficiency and lower mass. 
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Outpxit Characteristics of Thermionic Power Sources 
(a) Typical I-V Curve of Cylindrical Thermionic Con 
verter (h) Linear Approximation Representing the 
Output of a Thermionic Reactor. I = 2U390 amps. 

V = 77.8 volts ^ 






Fig. h Self-Field Magnet oplasmadynamic Acreleretor (MPD Thruster) 
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The MPD thruster can use a variety of propellants since any medium 
through which an arc can be passed is capable of being accelerated. The pro- 
pellant choice depends on the mission. Typical experimental MPD thrusters, 
currently being tested, use inert gases as propellants. 

The voltage- current characteristic of an experimental MPD thruster 
previously tested at Princeton (Ref. 5) is shown in Fig. 5. As shown in the 
figure, there are two branches to i;he V-I characteristic. In the upper branch 
(high slope) uhe voltage is proportional to In order to model the V-I 

oiiaracterlstic of the thruster for analytical purposes, the following equation 
was used. 


= Vp + gl^ (3) 

In Eq. (3) the thruster voltage is proportional to the 2nd power 
of the current. The form of Eq. (3) was selected even thou^ a cubic approxi- 
mation would be better because a cubic V-I characteristic leads to mathematical 

complexity in subsequent equations. The analytical approximation with V = 30 

-6 ^ 
and 3 = 1*2 x 10 is compared with an experimental characteristic in Fig. 6. 

Ij.O INDUCTIVE COUPLING CIRCUITS 
U.l Transformer Type Circuit 

To provide maximum flexibility in voltage transformation, it might 
be desirable to have an inductive unit with separate primary and secondary 
windings, i.e. a transformer. A schematic diagram which shows a transformer 
type circuit for coupling a thermionic reactor with a MPD thruster is presented 
in Fig. 7. By proper choice of the turns ratio between the primary and secon- 
dary, any desired voltage at the output can be obtained. The initial phase of 
this study centered around this type circuit. 

Analysis shows that good coupling is possible at a variety of turns 
ratios. An important requirem^ ", for good efficiency is the speed of the con- 
trol cl^-r.-:ui . 3^ (see Fig. 7)- The control element must be capable of inter- 
rupting the primary current in a time significantly less than 1 msec. If 
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(volts) 



THRUSTER CURRENT (KA) 

Fig, 5 Performance of aji ExperimentaJL MPD 
Thruster at Princeton University. 
Propellant is argon. The depend- 
ence has a theoretical basis and does 
not depend on the electrode geometry. 


THRUSTER VOLTAGE (VOLTS) 



Fig, 6 Comparison of Analytical Approximation 
and Experimental MPD Tliruster Charac- 
teristic. The approximate equation is 
= Vp, + 01^, with Vp “ 50 volts and' 

3 = 1.2 X 10“^ (amps)“^. 
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Fig. T Transformer Type of Circuit for Coupling Thermionic 
Reactor and MPD Thruster 
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does not reach a high impedance state ("open circuit") rapidly, large amoxants 
of power are dissipated in the switch. 

1*,2 Self Inductor Circuit 

One analytical restilt obtained with the transformer circuit is parti- 
cularly interesting. This shows that a turns ratio of 1, i.e. a self inductor, 
is compatible with the nominal values stated earlier for the current and voltage 
in the thermionic converter array. A self inductor type of c\irrent is preferred 
because of the additional simplicity and lower mass available with a single 
winding. Consequently, greater emphasis has been placed on the self inductor 
circuit in the recent phase of this work. 

A schematic of the self inductor circuit is shown in Fig. 8. As shown 
in the figure, the thruster is connected across the switching unit in this type 
of circuit. The switch is closed during the charging cycle and the current 
by-passes the thruster to build up a magnetic field in the inductor. During the 
thrust part of the cycle, the switch is open and the current is then forced 
to flow through the thruster unit. The current decays rapidly during this time 
consoming the energy which had been stored in the inductor during the charging 
cycle. 

The particular device to be used as the switch in the circuit is not 
yet designed. However, an array of transistors or SCR's could be used for this 
purpose. For performance calculations it was assumed that a one volt drop 
occurs across the switch when it is in the closed position. A conceptual de- 
sign for a switch of this type could consist of 100-200 semiconductors in parallel. 
Each semiconductor carries ^ 100 amps. Previous system analyses for spacecraft 
power conditioning have used transistors for inverter switching with similar 
operating characteristics (Ref. 6). 

5.0 PERFORMANCE OF THE SELF-INDUCTOR CIRCUIT 
5.1 Current Buildup 

Using the linear approximation for the thermionic generating unit, a 
differential equation for the current buildup during the charging cycle can be 
written as follows : 
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Fig. 8 Self Inductive Ty]pe of Circuit for Coupling a Thermionic 
Reactor and a ^^PD Thruster 
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0 


ik) 


L 


dt 


+ I R 

P 



or 


^ + i Zo 

dt T ■ L 


0 


(5) 


where 


T 


R 

P 


V /I 
o o 


( 6 ) 


So-lvirtg- Eq. (5)» and evaluating the constants using initial condi- 
tions yields the following result. 


I 


I 



V 



o 


1 - exp(- 

R + V /I 


1 

o 

o 




(7) 


5.2 Current During the Thrust Pulses 

The thrust pulse hegins when the switch is opened. The switch remains 
open during the thrust pulse. The turn-off time of the switch is assumed to 
he negligible (« 1 msec). It may be possible to initiate the arc in the MPD 
thruster by the voltage transient as the switch opens. If this proves to be 
unfeasible, it is possible to use a small auxiliary pulse to start the arc in 
the MPD thruster. At the time the switch is opened, the current is at its maxi- 
mum value. The decay of the current during the thrust pvilse is described by 
the following differential equation. 


^ i * ''f * ei" - - f > “ ° 

o 


( 8 ) 
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During the thrust pulse the current falls to its final value I. 
Solution of Eq,. (8) yields a relationship for I as follows. 


T = T(..l - K exp(-i/q t*) 

I - I(t») - t») 


where 


_ 2c $ + b — i/^ 
2c 5 + b + »^ 


- .w2 1 

q = b - 4 ac 


(9) 


( 10 ) 


( 11 ) 


a = 



( 12 ) 


b 


R 


R 


+ R + V /I 
s o o 

L 


(13) 


ilh) 


The values of I, I, t* and T as shown in Eqs. (7) and (9) are sub- 
ject to optimization. It is desirable to have the current during the charging 
cycle remain as close as possible to the maximum power point of the thermionic 
converter array (~ 12000 amps). On the other hand effective use of the inductor 
requires that the difference between I and I be as large as possible. It is also 
desirable to maintain the current through the thruster as hi^ as possible. Be- 
cause of these considerations, it is necessary to evaluate I, I, t* and T as 
part of the comple+e optimization. 
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6.0 CHARACTERISTICS OF THE INDUCTOR 


6.1 Inductor Configuration 

A system optimization requires detailed relationships for the variables 
which affect the characteristics of the inductor. A toroidal type of inductor 
may be preferred because it completely contains the magnetic field. Time-varying 
stray magnetic fields may be a source of interference for spacecraft instru- 
mentation. However, it is difficult to design a toroidal inductor which will 
fit inside the shadow cone as shown in Fig. 1 . As a result, a solenoidal inductor 
geometry was chosen for this study. A solenoid also has a better ratio of induc- 
tance to resistance than a toroid of equivalent mass. The size of the inductor 
shown in Fig. 1 is approximately to scale for a near optimum solenoidal confi- 
guration. 

6.2 Inductance 

The dimensions of the solenoidal inductor are shown in Fig. 9- All 
dimensions are in centimeters. 
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The inductance of a solenoid of this type is given by Eq. ( 15 ) 

below: 

L(yh) = ^ jpd ^ (.693 B^)j (15) 

The functions F and B are tabulated in Ref, 7 . However, the following analy- 

8 

tical expressions have been obtained for these functions by curve fitting to 
the tabulated data. 


F = .02378 (d/h) - .00652 (d/h)^ 


(16) 


B 

s 


.28 - .44 e:q) 


I"- .45086 (h/s) 


( 17 ) 


6.3 Mass of the Inductor 

The mass of the inductor varies inversely with the resistance of 
the winding. It’s desirable to have the lowest possible resistance so it is 
necessary to con5>romise between the mass and the ohmic loss. Ihe material 
which has the lowest ratio of mass density to electrical conductivity is 
aluminum (sodium has a lower value but it is not a practical noice). Con- 
sequently, the inductor windings are made of aluminum which is laminated to 
cut down eddj’" current losses. The laminations are assumed to take up lOjJ of 
the volume. 

6.4 Resistance 

The resistivity of aluminxim is temperature dependent. Since signi- 
ficant ohmic losses will cause a temperature rise in the inductor it is neces- 
sary to use temperature dependent resistivity. An equation was fit to published 
data for the resistivity of aluminum (Ref. 8 ) as shown below, 

- cm) = 5.915 + .0232 e - 1.75 X 10”^ 0 ^ (18) 
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The value of the resistivity given by Eq. (l8) is increased by lOjC to account 
for the volume lost by lamination. This is then used to calculate the resis- 
tance of the helical inductance winding in the following equation 


R 


M (TTd)^ (N> 1)^ 
- h^ 



(19) 


6.5 Heat Balance 

The thermal power generated by ohmic losses in the inductor is radiated 
from its outer surface. No additional radiator is necessary. The temperatvire 
of the inductor rises to approximately 600 % at thermal equilibrium by radia- 
tion. The above relationships for the characteristics of the inductor combined 
with previously derived equations for the performance of the thruster and the 
thermionic generating unit contain the necessary elements to perform a system 
optimization. 


7.0 OPTIMIZATION 

7. 1 Coupling Effectiveness 

Two criteria were used for optimization of the system. One is to 
maximize the coupling effectiveness, the other is to minimize the total spe- 
cific mass of the power generating system and thrust unit. 

The notion of coupling effectiveness requires some definition and 
discussion. The coupling effectiveness is defined as the average net power 
reaching the thruster divided by the maximum power available from the thermionic 
generating unit. This can be written as: 


n 





Energy Per Thrust Pulse 
T + t* 


') 


( 20 ) 


This is not precisely the same as an efficiency. A major reason 
why n is not unity is that the ciurent varies during a cycle and does not always 
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stay at the maximum power value, 1^/2. Only a portion of the power which is 
not coupled into the thruster is dissipated in ohmic and switching losses. 

For example, neai the optimum approximately lk% of the maximum available power 
is not transferred to the thruster (n opt cs %3%)\ ’’•-ut only approximately 6j{ 
is consumed by ohmic losses. 

Thus, the value of ri is a correct measure of the generating capacity 
required by the system, but it does not correctly determine reactor thermal 
power because the thermal input to the thermionic converters varies with current. 
The specific mass of the reactor is inversely proportional to the coupling ffec- 
tiveness, but the specific mass of the radiator does not vary in the same way 
and neither does the fuel burn-up. Because of these facts, the estimates of 
specific mass for the entire system produced in this study are eomewhat con- 
servative. 


7.2 Specific Mass 

A previous study (Ref. U) estimated the specific mass for the therm- 
ionic heat pipe reactor as a = 19.7 kg/KWe. This value was based on a system 
where the reactor gross power is 1»7^ KWe and the net power is UOO KWe. It's 
convenient to redefine the specific mass value in terms of the gross power, 
thus 





l6.6 kg/KWe 


( 21 ) 


This specific mass of the switching device is difficult to estimate 
since the switch is not defined. However for calculation purposes, the mass 
of transistors previously used for spacecraft power conditioning (Ref. 6) will 
be used. Each transistor caries 75 amperes; hence the number of transistors 
required is 


Number of Transistors 


75 



75 

o 


2 P 

75 V^n 


( 22 ) 
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Ihe mass associated with the transistors and the appropriate connections is 
■ 35 kg per transistor (Ref. 6). Multiplying by the number of transistors given 
by Eq. (22) gives the specific mass of the switching device, = .12 kg/KWe. 

A small auxiliary radiator is required to dissipate the heat generated 

in the switching transistors. Specific mass of this radiator, a , is estimated 

In 

to be about 1 kg/KWe (Ref. 6). 

The specific mass of a quasi-steady MPD thrusLer has been estimated 
as = .6 Kg/KWe in the power range of a few megawatts (Ref. 9)» The total 
mass of the system can be written in terms of the net power and the coupling 
effectiveness as follows 


a_ P a_ P a__ P 

R T . TR . 

+ — + 


n 


n 




P - Mass of System 


(23) 


Dividing through by the net power, P, yields the specific mass of the entire 
system 


a . 


OL, 

^ T TR 


+ a, 


M 



i2h) 


or 


a 


S 


17. T ^ 
n 



(25) 


Inspection of Eq. (25) shows that a varies inversely with the 
coupling effectiveness n and directly with the mass of the inductor The 

variation of the thruster efficiency is not directly included in the optimization. 
Hence, an arbitrary decision was made to keep the thruster current above 9j000 
amperes for efficient operation. A nominal 1 msec duration was chosen for the 
thrust pulse which is near the optimum. Other variables were then optimized 
to yield a minimxjm in the system specific mass by an iterative calculation using 
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a HP 982 OA prograinmable calculator. This calculator program and sample output 
are included as Appendices A and B respectively. 


8.0 OPTIMIZATION RESULTS 

The optimum value fox the parameters of the inductive coupled system are 
as follows : 

• Duration of the thrust pulse, t* = 1.0 msec 

e Time between pulses, T = lf.5 msec 

• Maximum current I = 1^,7 kA 

• Minim’mi current I = 9*1 kA 

• Inductance, L = 29 yh (9 turns) 

• Mass of the inductor M^ = 950 kg 

• Coupling effectiveness n = 86$ 

• S^^stem specific mass ot = 23.5 kg/KWe 

O 

Variations of the parameters in the vicinity of the optimum are shown in 
Fig. 10, 11, 12 and 13. The design point for the inductor mass (950 kg) is 
slightly higher than the value which yields a minimum in the specific mass. 

This value was selected to yield a higher coupling effectiveness and also to 
provide some safety factor in the design. The performance has a broad optimum 
with inductance in the range of 15-60 Vh (7 to 13 turns) as shovn in Fig. 11. 

The optimum value of the pulse duration is slightly above 1 msec as shown in 
Fig. 12. However, the 1 msec value was chosen as a design point to maintain 
a higher value of I (refer to Eq. (9)). Similarly the optimxmi value of the 
time between pvilses, T, is slightly less than U,5 msec; however, the U.5 msec 
value was retained as a design point to give the value cf I (refer to Eq, (?))• 

The waveforms of the cxirrent through the inductor and the power produced 
by the thruster are shown in Figs. ll+(a) and (b). The duty cycle of the thruster 
is 22$ at the design point. 
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Fig. 11 Variation of Specific Mass and Coupling Effectiveness With 
Inductances, L. Other parameters not varied in this figure: 
t* = 1.0 msec; T = lt.5 msec; = 950 kg. 
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Fig. 12 Variation of Specific Mass and Coupling Effectiveness With 
Pulse Duration, t*. Other parameters not varied in this 
figure: T = ii.5 msec; L = 29 yh; M_ = 950 kg. 
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Fig. 13 Variation Specific Mass and Coupling Effectiveness With 
Time Between Pulses, T. Other parameters no"; varied in 
this figxxre: t* = 1.0 msec; L = 29 yh, = 950 kg. 
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Fig. lU Waveforms of Inductor Curj^ent and MPD Iliruuter Power at 

Optimized Design Point. I = lU.T kA; I = 9-1 kA; T = msec; 
t* = 1.0 msec. 
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9.0 SUMMARY AND CONCLUSIONS 

The results of the preliminary analysis have shown that an inductively 
coupled thermionic reactor and an MPD thruster is a potentially attractive 
system for spacecraft electric propxilsion. Good coupling efficiency has "been 
shown to he possible. Specific ma.ss of a nuclear electric propulsion system 
using ion rockets is calculated to be ~ 28 kg/KWe (Ref. 3). Thus the specific 
mass estimate of 23-5 kg/KWe for the TI-MPD system obtained in this work appears 
to be competitive. 

Additional experimental and analytical work must be accomplished before 
the system is proven flight ready. Experimental data are needed to assess the 
performance of an MPD thruster when driven by an inductive source. An important 
remaining design uncertainty is the current interrupting switch (solid state or 
other). Design of a switching circuit to interrupt large currents in an in- 
ductive circuit is somewhat difficult. 

If the remaining considerations can be satisfactorily treated by continued 
development, the thermionic - MPD system will be an important option for future 
space missions. 
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APPENDIX A 


A COMPUTER PROGRAM FOR CALCULATING 
THE PERFORMANCE OF THE TI-MPD SYSTEM 
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APPENDIX B 


TYPICAL COMPUTER OUTPUT - NEAR OPTIMUM CONDITIONS 
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